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Abstract 

We perform a spectral decomposition of the dynamical correlation function of the 
spin 1/2 XXZ model into an infinite sum of products of form factors. Beneath the four- 
particle threshold in momentum space the only non-zero contributions to this sum are 
the two-particle term and the trivial vacuum term. We calculate the two-particle term 
by making use of the integral expressions for form factors provided recently by the 
Kyoto school. We evaluate the necessary integrals by expanding to twelfth order in 
q. We show plots of S(w,k), for k = and n at various values of the anisotropy 
parameter, and for fixed anisotropy at various k around and n. 
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1 Introduction 

In the last two years, a remarkable series of papers have appeared in which the consequences 
of the quantum affine symmetry of quantum spin chains is explored [||, |2|, |3], ||, ||. In the 
simplest case, the spin chain is the antiferromagnetic spin-1/2 XXZ model, and the associated 
non-abelian symmetry that of U q (sl(2)) [[l], |^]. The U q (sl(2)) symmetry manifests itself in 
different ways. Firstly, the Hamiltonian of this model commutes with the algebra directly. 
Secondly, the model has a vertex operator algebra as a dynamical symmetry. These vertex 
operators are in turn intertwiners of certain U q (sl(2)) modules. However, we do not wish to 
review the work of the Kyoto School here. An excellent and accessible reference is the recent 
book by Jimbo and Miwa |J. It is sufficient for our purposes to note that this approach 
leads to exact expressions for both correlations functions and form factors of local operators 
of the XXZ model. These objects are given as the trace of vertex operators over irreducible 
highest weight representations of U q (sl(2)). Using the free field representation of the algebra 
these traces may be evaluated to yield integral expressions. 

Perhaps somewhat surprisingly, one-dimensional quantum spin chains are relevant to 
certain physical systems (crystals that have a much stronger interaction in one direction, 
making them effectively one-dimensional - for reviews of theoretical approaches and exper- 
imental data, see ]7], [|| and references therein). A quantity of key interest in this context 
is the dynamical correlation function (defined in the next section). This quantity is related 
to neutron scattering data |J. In this paper we calculate the dynamical correlation of the 
spin-1/2 XXZ model by performing a spectral decomposition, i.e., we re-express the corre- 
lation function as an infinite sum over products of form factors. Since we are calculating 
the dynamical correlation function in momentum space, all terms other than the vacuum 
and two-particles one vanish if we restrict ourselves to beneath the four-particle threshold 
(odd particle terms are always zero). This may seem an artificial restriction, but we are 
lucky in that experimentalist are, for the most part, only interested in the region around 
the two-particle threshold. The vacuum term is trivial to evaluate. It gives a delta function 
peak at w = and k = ir. 

We evaluate the two-particle form factor by expanding the integral expression of the 



Kyoto School in q (to twelfth order). Up to some subtleties about the integration contour 
(which we explain), the coefficients are then relatively easy to work out. 

In Section 2, we define the XXZ model, and describe the spectral decomposition of the 
dynamical correlation function. In Section 3, we write down the form factor expressions 
derived by the Kyoto School. We describe the general method of q-expanding such integrals, 
and evaluate the specific integrals required to twelfth order in q. In Section 4, we present the 
results of some consistency checks on our form factors. We calculate all the contributions 
to the correlation function, and show graphically how the correlation function behaves in 
various regions of parameter space defined by energy, momentum and q. Finally, in Section 
5, we draw some conclusions. 

2 The Dynamical Correlation Function of the XXZ 
Model 

2.1 The XXZ Hamiltonian 



The Hamiltonian of the spin 1/2 XXZ Heisenberg quantum spin chain is Jl0|, [II], |12 , 



1 oo 

W = -o £ (°*°i-i + °*°»i + A(T ^+i)> (2- 1 ) 

i=— oo 

where A = (q + g _1 )/2 is an isotropy parameter. Here o{ are the Pauli matrices acting at 
the i th site, with the Hamiltonian acting formally on the infinite tensor product 

W = ---V®V®V®V---, (2.2) 

where V ~ C 2 . We consider the model in the massive antiferromagnetic phase A < — 1, 
which we parametrise by — 1 < q < 0. 

2.2 The dynamical correlation function 

The dynamical correlation function we consider is 

oo 

S*(w,k) = J dtE e*( v *- k P*S%p), 

-oo P ez (2.3) 

S l (t,p) = i< vac\a(t,p) • a(0,0)\vac>i, 



where a^(t,p) refers to <ji acting on the site at position p at time t (in Minkowski space), 
and i = 0, 1 refers to a choice of one of the two antiferromagnetic boundary conditions. 
We may define these boundary conditions by arbitrarily choosing the eigenvalue of a 3 at 
the j'th site, s J = (— 1) J+J , for \j\ >> 1. The dynamical correlation function is of physical 
interest because of its connection with neutron scattering data || || . For a physical system 
described by the Hamiltonian |2.1| , the cross section for scattering of neutrons that give up 
energy w and momentum k is proportional to S t (w, k). In this paper we describe a method 
of calculating this function as an expansion in the parameter q. 

2.3 Decomposition of the identity 

We evaluate |2l| by performing a spectral decomposition (a similar method for calculating 
the dynamical correlation function of the s = 1 XXX model, as approximated by an 0(3) 
non- linear sigma model, was used in reference [HJ). That is, we insert a complete set of 
energy eigenstates between the Pauli matrices. The Kyoto School has shown how to construct 
these states within the context of the representation theory of the quantum affine algebra 
U q (sl(2)). We don't wish to go into excessive detail about this rather long story, but rather 
pull results from the literature as necessary. A complete discussion of this subject is given 
in reference ||. 

Jimbo and Miwa have conjectured (and given strong plausibility arguments for) the 
completeness of the following decomposition of the identity in terms of the eigenstates of the 
Hamiltonian |J: 

I=EE E — id£l---d£n\£m---i£\ >e n ,-,e 1 ;i i;e 1 ,-,e n < ^1, ' ' ' An\, (2.4) 

i=0,ln>0e n ,-,e 1 n - J 

where, d£ = d£,/(2niz), the contours are around the unit circle, |£j| = 1, and the pseudopar- 
ticle 'charges' are Cj = ±.1. 

The action of the spatial shift operator and Hamiltonian on multiparticle eigenstates 
states is given by 

n 

*t (2.5) 

Hxxz\€l, " ' ')£n >i = E -E(£i)l£l) ' ' ' i£n >i, 



where, 

r(0 = i 



-i y(q£ 2 ) 

OAt-'V (26) 
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m = ^£iogT(o, 



2g s d£ 



DC 



and 6> a (6) = (a; a)oo(6; a)oo(afe 1 ;a) 00 ; (a; 6)oo = II (1 — a& n ). For later purposes, we also 

n=0 

oo 

define (a; 6, c) = II (1 — ab n c m ). As before, z = 0, 1 corresponds to a choice of one of 

n,m=0 

the boundary conditions. The states |2.5| correspond to the hole states in the Bethe Ansatz 
spectrum 



Inserting this expression for the indentity into |2.3| , and making use of the properties, 



a j {t,p) = e mf T-?V'(0,0)TPe- iWi , 
T\vac >i = \vac >i_ 



-i ■ 



we obtain the following expression for the correlation function: 

S l (w,k) = J2 S l n {w,k), where, 

n>0 

S n (w, k) = E JC jft h. § *i • • ■ din ( T{ ^ n) ) P S(w - Efa) ■■■- E(z n )) (2.8) 

x i+p < vac\a(0,0)\£ n , ■■■,£,! > £ni ... i£i;i+ p i ;ei ,..., en < &, • • • , fn|£(0, 0)\vac>i, 

where the boundary conditions are understood as modulo 2. Using the delta function A (z) = 
E z p , we may rewrite this expression as 



S l n (w,k)= E A ^^1 • • • ^» A ( r(ft)Ctt„)0 *( w ~ S(fr) • • • ~ £(*»)) 

X (i < wac|a(0, 0)|f„, • • • , 6 >e n ,-,e 1 -,i i-, ei ,-,e n < £i, - - - , fn|£(0, 0)|vac >i 

+ ( r(6) e - fc r(€n) ) I"* < WaC k( ' )Kn' ' " ' > 6 >e n ,-,«i;l-i i;ei,-,£„ < £l, ' " ' , fn|oi(0, 0)\vdC > t 

(2.9) 
Two of the integrals may be carried out by making use of the identities, 

$ (£%) mOMO = g^flt> ttol/Ko) = 0, |£o| = 1}, 

K4)a(/(omo = g ^ffiU, ' {wee,,) = i, i&i = i}. 

The first observation we make about |2.9| is that, for a given choice of w and k, the 
integral is only non-zero if their exist £1, • • • , £ n such that both, t(£i) 2 • • • i~(£ n ) 2 = e 2lk , and, 
-E'(^i) + • • • + E(£ n ) = w hold (that is, there is momentum and energy conservation). In 
particle physics language, w needs to be greater that the n-particle threshold. (If we were 



dealing with free relativistic bosons, the analogous argument would imply that the energy 
and momentum conservation conditions could be met only with w 2 > k 2 + n 2 m 2 .) Putting 
the argument another way, for a given choice of w and k, the series [2.3| terminates. 

Charge conservation in this model implies that S^ odd (iu, A;) vanishes. (This point is 
clear within the context of the Kyoto School approach when form factors are expressed as 
traces over vertex operators. Ultimately, the charge conservation is related to the charge 
conservation in the Boltzmann weights of the corresponding six vertex model ||15|| .) So up to 
the four-particle threshold, the is an exact equality between S l (w, k) and Sq(w, k) + S l 2 {w, k). 
The existence and location of the two and four particle thresholds are shown in Figs. 1 and 
2. In Fig. 1, we simply choose a large number of (£1,^2) pairs randomly destributed over 
|£i| = 1, and plot w = -E'(^i) + Efa) versus k = — ilog(r(£i) +r(£ 2 )) (we evaluate E(£i) and 
r te) by the q expansions given in Section 4). In Fig. 2, we do the same but with (£1, ■ • ■ , £4). 
The thresholds are the lower edges of the two scatter plots. 

The vacuum, or zero-particle, term, is given by 

S l (w, k) = A(e 2ik )5(w) L < vac\a(0, 0)\vac > t +e ifc i-i < vac\a(0, 0)\vac >i_ij 

(2.11) 

i < vac\a(0,0)\vac>i . 
Again, charge conservation implies that j< vac\a ± \vac >i= 0. Using a Z 2 symmetry given 
in §; 

i < vac\a z (0, 0)\vac >i= - i_j < vac\cr z (0, 0)\vac >%-i, ( 2 -12) 

we can simplify this expression to give, 

S l {w, k) = A{e 2ik )5{w) G < vac\a z {0, 0)\vac >*) 2 (1 - e ik ). (2.13) 

So the vacuum contribution to S*(w,k) is a delta function peak at k = it. The co- 
efficient i < vac\a z (0, 0)\vac >« is non-zero, and is in fact the staggered polarization = 
(# 2 ; Q 2 )lo/(~ < l 2 i Q 2 )lo derived by Baxter [0, and reproduced by the Kyoto School (by explic- 



itly integrating the relevant integral formula pj). 



In order to calculate the two particle contribution, we use the identities |2.10| to carry out 
the two integrals. This gives, 

£2,ei 6,6 Ul,42j (2.14) 

+ ( r(giM$„) ) 1-* < vac h\&Z,l >6 2 ,ei;l-i i;ei,ea < 6> 6k| uac >i 



Here, the sum is over all (6, 6) which are solutions of t(6) 2 t(6) 2 = e 2lk and E(£i)+E(£ 2 ) = 
w. The Jacobian factor c(6,6) is given by 

c(£i,£ 2 ) = -^|6£'(6)£(6) -6£'(6)£(6)l- (2.15) 

g-g 



3 Evaluation of Form Factors 

In this section we shall give expressions for the form factors appearing in equation [2.14 . 



3.1 The results of the Kyoto School 

The Kyoto School have shown how to calculate form factors of the XXZ model as traces 
of vertex operators over infinite-dimensional irreducible highest weight modules of U q (sl(2)) 
- see |J and references therein. By bosonizing the vertex operators, and representing the 
highest weight modules in terms of a bosonic Fock space, they are able to write down explicit 
expressions (generally given in terms of multiple integrals of elliptic functions) for both form 
factors and correlation functions. The following expressions for the form factors of interest 
are given in reference ||: 

t <vac\a |6,6 >-,-,= (-?) 6 6 (g W ) 4 p(g) A - q 3 /ui)9A - q z /u2) > ( 3 ' 16 ) 

4 <^c|^|6,6 >+,-,= (-g)"*(i - g- 2 )(g 2 )!(g 4 )|p(g) 2 7(« 2 M) 

x § dv § dw(qv~ l — q~ l Ui l )w l ~ % (—v) % (—q~ l v)2(—q ? '/v)2 fS 17") 

C C++C- \ ■ J 

x q 8(-q 4i v 2 /(w 2 U1 U2)) e k - l (-qUkM4(-q 3 w/u k ) 

(w)2(q 2 /w)2(-q- 1 v/w) 2 (-q 3 w/v)2 k ±* 2 (-qu k ) 4 (-q 3 u- 1 )i{u k /v) i {q~ 2 v/u k ) i ' 

where Ui = — 6 2 - We introduce the notation (a) n = (a; g n )oo and (a) nim = (a; g n , q m )oo- 7 
and p are defined by 



7(a) 

p(g) 



(3.18) 



(q 6 u)4 t 4,(q 2 u 1 )4,4' 



(9 6 ) 



14,4 

.6^ 



The contours C, C± are denned by 



for n > 0, 

for n > 0, and, 

for n > 0, or equiv., 

-wq 2n+i < v < —wq~ 2n+l for n > 0. 



C: 


Mjg 4ri 


< U < 


Ui q 2 - 4n 


C±: 


g(2n+l)±l 


< W < 


9 -(2n+l)±l 




2m— 1 

—vq 


<W< 


_ v q-2n-3 




„,.„2n+3 


^ ii <^ 


„,.„-2n+l 



(3.19) 



We use the >, < signs to imply that the contour on the left-hand side runs outside or inside 
a series of poles. The contour C is not a circle, because \uiq 2 \ < |wj|, but it is nevertheless 
well defined - it simply loops inside the pole at u,iq 2 . 

The remaining non-zero form factors are defined through the Z2 symmetry of the theory 
I by 

i<vac\a-\£ 2 ,£i>++;i = i_i< vac\cr + \Z,2,£i >— ; i-», . . 

i<vac\a z \£ 2 ,£i >-+;i = -i-i< uac|o- z |^ 2 ,6 >+-;!-« • 



All non-zero dual form factors are also defined in terms of |3.16| and |3.17| via the symmetry 



i; eil e 2 <£l, 6|CkoC >i= i< vac\0\ - q£ X , -q& >-ei,e 2 ;i, ( 3 - 21 ) 

where represent any local operator in the theory. 

3.2 Evaluation of the integral 

To obtain a useful expression for the dynamical correlation function |2.14|, we need to evaluate 



the integral |3.17| . It is a daunting task to attempt this directly, and we instead employ the 



technique of q expanding the integral. In order to illustrate some of the subtleties of this 
technique, we apply it first to a simpler integral I(k, I), defined by 

JM= /*(K5bo? (3 - 22) 

where the k and / are integers, and the contour is chosen as C : q l+An < £ < q~ k ~ 4n f or 
n > 0. We consider three generic cases. 



k,l > 

In this case we can simply expand each term in the products in the denominator as a series, 



i.e., 



I(k,l) =/^IIE (£g fe + 4 ") mi £ (rV+ 4n r 2 - (3-23) 

f~< n>0 mi>0 m2>0 

We may rewrite this as, 

I(kJ)=id^an(Oq n , (3-24) 



£ n>0 

where for any n, a n (£) is a finite polynomial in £. We calculate the residue at £ = by 
extracting the constant term in a n (£). As an example we find, 

1(2, 2) = 1 + q A + 3q 12 + 12g 16 + 21q 20 + 38g 24 + 63g 28 + 106g 32 ■ • • . (3.25) 

k< 0,1 > 

Consider the example J(— 1,3), where C : q 3 < C, < q. Now, we cannot expand the denomi- 
nator in positive powers of g as \^/q\ < in the (1 — £/g) term. To be able to do so we must 
shift the contour outside of the £ = q pole, i.e., 

/( - 1,3) = "/% - «/o(W/o. + rai' (3 ' 26) 

where C : q < £ < q~ 3 , and the second term comes from the residue at £ = q. We can then 
calculate the integral over the contour C by the method used in the previous case. 

k> 0,1 = 0, or k = 0,l > 

In each of these cases we can q expand the integrand without problem, but we also obtain a 
pole at £ = 1 which we must exclude or include for the two cases respectively. 



We now show how we have applied this q expansion technique to the integral |3.17| . The 
location of poles in the integrand is determined by the following product of terms in the 
denominator: 

{w) 2 {q 2 /w) 2 (-q~ l v/w)2{-q*w/v)2{ui/v) A (u2/v) i (q 2 v/ui) i (q~ 2 v/u2)A. (3.27) 

We fix w at a point in the C + or C_ band (C± are genuinely bands, unlike C) and carry 
out the v integration first. In order to q expand we must shift the v contour outside of the 



v = q 2 U2 and v = —qw poles, in order that it lies within the bands \ui\ < \v\ < \q~ 2 Ui\ and 
\qw\ < \v\ < \q~ 1 w\ (we call this contour C). The entire integrand may then be expanded. 
In the v integration of this expanded integrand we are left (at any order in q) with poles at 
v = Mi, u 2 , 0. In the subsequent w integration we get poles at w = 0, and also at w = 1 for 
the C- contour. It remains to calculate the w integrals associated with the residues at the 
v = q 2 U2 and v = —qw poles. At v = q 2 u 2 we get the following w poles in the denominator, 

{w)2{q 2 /w)2{-qu 2 /w) 2 {-qw/u2)2- (3.28) 

However the third factor is partially cancelled by a term in the numerator leaving us with, 

{w) 2 {q 2 /w) 2 {-q*u 2 /w) A {-qw/u 2 ) 2 . (3.29) 

In determing the location of the w contour with respect to the U2 dependent poles we 
encounter a new subtlety. The origin of the rather curious original C contour was that 
it arose from the sum of two expressions with contours C\ : \q A Ui\ < \v\ < \q 2 Ui\ and 
C-2 : | itj| < \v\ < \q~ 2 Ui\ [[J. (We don't wish to go into the origin or form of these two 
expressions - yet again we refer the reader to reference 0.) Compatibility of these contours 
with \wq 3 \ < \v\ < \wq\, requires \q 4 Ui\ < \wq\ < \ui\ and \ibi\ < \wq\ < \q~ 4 Ui\ respectively 
- which are of course incompatible with each other. So in order to work out whether we 



should choose \w\ < \u 2 /q\ or |w| > \u 2 /q\ in the integral associated with |3.29| , we must 
determine whether the v = q 2 U2 arises from the part of the original integral associated with 
the contour C\ or the part with contour C*2. In fact it arises from the C\ part and so we 
must choose \wq\ < \v,2\- This is convenient, because it means that we can carry out the q 
expansion without shifting the contour any further. 

Now consider the w integral associated with the residue at the v = —qw pole. In this 
case we are left with the following w dependent terms in the denominator, 

(-u 1 /(qw)) 4 (-u 2 /(qw)) 4 . (3.30) 

Again, other products cancel with those in the denominator. We find that these poles arise 
only from terms associated with the C2 contribution to the original integral. This means 
that we must choose \ui\ < \wq\. This causes a problem because we cannot then expand the 

9 



(1 + u 2 /(qw)) an d (1 + ^/(qw)) factors in positive powers of q. Hence we must shift our 
contour to \qw\ < \ui\, which requires that we evaluate the residues at two extra poles. 

A complete list of the combinations of poles which we must evaluate in order to obtain a 
complete q expansion for the form factor |3.17| is given in Table 1. The sign (±) in the table 
indicates the sign of the contribution of a particular residue to the overall integral. 



Table 1 - Location of residues evaluated in association with the form factor 3.17 



V 


w 





(+)0 (+)1 


1 


(+)o (+)i 


u x 


(+)o (+)i 


u 2 


(+)o (+)i 


q 2 u 2 


(-)o (-)i 


—qw 


(-)O {-)-u x /q (-)-u 2 /q 



The technique for calculating i<vac\a z \ — q£, —q£ >+_;, is completely analogous (although 
not identical as different poles contribute). 

4 Results 



In practice, as a non-trivial check of our technique, we calculated each of the following eight 
integral form factors separately: i<vac\a z \£, 2 , £i >+-;*> i< vac \°' z \^2,(,i >-+;i> 
i<vac\a z \ — q£i, —q£,2 >+-;« an d i<vac\a z \ — q£i,—q£, 2 >-+-,% (each for i = 0,1). We did 
this to 6th order in q, and found that the Z 2 symmetry relations |3.20| held. In addition we 
found that the dual form factors were precisely the complex conjugates of the form factors, 
i.e., 

i<v ac\a z \ - g£i, -g£ 2 >- ei -e 2 ;i= i<'uac|cr 2r |^ 1 , ^ _1 >e 2 ,er,i ■ (4.31) 



This equality is also true for the exact expression 3.16, i.e. 



i<vac\a + \ - g£i, -q£ 2 >_ 6li 



-e%i t 



<vac\a \£ 2 \£ 1 >e 2 ,ei;i ■ 



(4.32) 



10 



Having established these identities (to sixth order in q at least), there are only two in- 
dependent sets of form factors (the members of each set being related by Z 2 symmetry 
or conjugation). We then proceeded to calculate one member from each set, in practice 
o< fac|<r 2 | — g£i, — g£ 2 >+- ; o an d i< fac|<r 2 | — q£i, — g£ 2 >+- ; i, to twelfth order in q. The 
results are given in Tables 2 and 3 at the end of the paper. 

Using these results, we calculated the two quantities required in the a z contribution to 
2TT4| (to order g 12 ), 



/f 7 (^i>6) = i< vac\a z \£ 2 ,ii >+-,i i< vac\a z \ -q€i,-q£ 2 >+-, 
+i< vac\a z \£ 2 ,£i >-+;% »< vac\a z \ -q£i,-q& >- 



(4.33) 
+i_i< vac\a z |f 2 ,fi >-+;i-i i< vac\a z \ -q£i,-q& >-+;* • 



/f JI (£i> 6) = i-i<vac\a z \£ t2 ,£i>+- ) i-ii<vac\cr z \-q£i,-qt;2> 



f zI and /* are both actually independent of i because of the Z 2 symmetry. This is also 
true for the exact a + and o~ terms, 

ft* fait*) = i<vac\a ± \&,£i> TT ;ii<vac\(T T \-q£ 1 ,-q&>±±;ii 
ft H (€u&) = i-i< vacl^ |£ 2 ,6 >TT;i-i i< vac\a T \ - <?£i, -g£ 2 >±±;t ■ 

Thus S 2 (w, k) is also independent of i. Furthermore we find, from our explicit results that, 

//"(£i,6) =/f 7 (±6,T6) = /?'&,&), ,,_ 

(4.35) 

/f 7/ (a,6) =-/f /7 (±6,T6) = #"(&,&). 

Again, the same properties hold for the corresponding a^ form factors. 



This later fact is important when we look more closely at formula |2.14| , which written in 
terms of the above functions, becomes, 

S l 2 (w, k) = c E 4^) ( 2 (tf '(fc.60 + /- z (6,6)) + /f J (£i,k)) 



+ Wfe) ( 2 (#■"&.&) + /*'"&, e 2 )) + #"(&,& 



It would be nice to able obtain S 2 (w, k) without performing the rather tedious task of solving 
the energy and momentum constraints 

r(a)M6) 2 = e- fc , 

£(6) + £(6) = ™, 
11 



to get all the £1 and £2 roots for a given w and k. This can be done, if given one pair 
of £1 and £2, that give a particular (w,k) (via the energy relation and t(£i)t(£i) = e~ lk ), 
one can get all the other solutions of the more general constraints [4.37| . From the ex- 



plicit form of r(£) and £?(£) (either |2.6| or |4.39| , we find that the set of such solutions is 



(£1, £2), (±£1, T-*^), (—£2, — £1), and another four with £1 <-> £ 2 . Given, the properties |4.35| , 
and the fact that c(£i,£2) = c(±£i,=f£ 2 ) = c(£2,£i), we can label S\{w,k) in terms of the 
original (£1,^2) pair, and write, 



53(6,6) = ^sy (2(/+ j (6,6)H-/r J (a,6)jH-/na,6) 

+ 2 (#"&,&) + /r /7 (a,6)) + /f 7/ (a,6)) . (4.38) 

For possible applications to neutron scattering data, it is useful, to work out Sl(w,k) 
as a function of w for fixed k values. We do this by numerically solving the momentum 
constraint t(£i)t(£ 2 ) = e, lk for a range of (£1,6) pairs, and calculating the associated w 
value and S % 2 (w,k) for each pair. In order to do this we identify r(£j) = e~ %p ^ di \ £j = e l6i , 
and use suitably truncated versions of the following g expansions of the functions |2]6] (we 
calculate c(£i,£ 2 ) i n the same fashion): 

P(*) = ^ + 2 E tSf 

m> r ° >! (4.39) 

In Figs 3 and 4, we plot S^iu, fc = 0) and S^w, k = n) for a range of g values. We find, 
in each case, that there is a rounded peak above the two-particle threshold. The location of 
the threshold moves to smaller w at q decreases towards — 1. This is as expected since, from 
4.39|, their is a mass gap in the theory, 



E( P = 0) = {q ~ q ~ l) V r Aqm „ (4.40) 

1 J 2 ^ m(l + g-)' l ' 

which vanishes as q —* — 1. The results at /c = and k = it are similar in form, although 

different in scale. The results at fc = and k = 11 are different since ^(±£1, =f£ 2 ) 7^ ^(^i, £2) 



- as follows from equation |4.38| . Our expressions for 5*2 (w, k = 0) and S^w, A; = it) are easily 
convergent at all three q values. As \q\ increases, the peaks become infinite spikes at, or very 
close to, the threshold (at least on the scale of the plots shown). 
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In Figs 5 and 6, we plot S2 (w, k) for a range of values at around k = and k — it, and at 
fixed q = —0.2. The location of the peaks move, as the location of the two-particle threshold 
moves as in Fig. 1. A rather curious fact is that the peak increases in height as k is increased 
from zero, and decreases in height at k is increased from n. 

5 Conclusions 

One of the criticisms sometimes made of the recent work on quantum spin chains, is that, 
despite all of the mathematical insights gained along the way, in the end one is left with 
formulae for correlation functions and form-factors which are intractable. We hope we have 
shown that, at the very least, these integral formulae can be regarded as generating expres- 
sions for q expansions of physically interesting quantities. We also hope that we have shown 
how, by making use of the expressions for form factors, one can use this approach to anal- 
yse dynamical quantities - rather that the short distance equal-time correlation functions to 
which it has formerly been applied. Of course, we have also calculated a quantity which, 
closer to q — — 1 may be relevant to the anisotropic spin-chains found in nature. We hope 
to comment on such applications in the future. 
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Figure Captions 

Fig. 1 A scatter plot of w = E(£i) + E(£ 2 ) vs k = — ?log(r(^i)r(^ 2 )) for (£1,^2) randomly 
distributed over < 6 t < 2n (where ^ = e %e% \ The lower limit of the points indicates the 
location of the two-particle threshold. Only the < k < 2n portion of the plot is shown. 

Fig. 2 A scatter plot of w = E(^)+E^ 2 )+E(^)+E(^) vs k = -1 logH&M&M&Mk)) 
for (Ci, ^2, ^3, ^4) randomly distributed over < 6i < 2n (where £j = e tei ). The lower limit of 
the points indicates the location of the four-particle threshold. Only the < k < 2-k portion 
of the plot is shown. 

Fig. 3 S 2 (w, k = 0) vs w for a range of q values. 

Fig. 4 S 2 (w, k = 7r) vs w for a range of q values. 

Fig. 5 S 2 (w, k) vs w for a range of k values close to k — 0, and for fixed q = —0.2. 

Fig. 6 S 2 (w, k) vs w for a range of k values close to k = n, and for fixed q = —0.2. 
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Fig. 2 Four Particle Threshold at q=-0.2 
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